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We study how the charge neutrahty affects the phase structure of three-flavor PNJL model. We 
point out that, within the conventional PNJL model at finite density the color neutrality is missing 
because the Wilson line serves as an external "colored" field coupled to dynamical quarks. In 
- - - this paper we heuristically assume that the model may still be applicable. To get color neutrality 

00 I one has then to allow non vanishing color chemical potentials. We study how the quark matter 

. phase diagram in (T, ma/^)-plane is affected by imposing neutrality and by including the Polyakov 

>— ^ ' loop dynamics. Although these two effects are correlated in a nonlinear way, the impact of the 

Polyakov loop turns out to be significant in the T direction, while imposing neutrality brings a 
5-H , remarkable effect in the m1/fi direction. In particular, we find a novel unlocking transition, when 

the temperature is increased, even in the chiral SU{3) limit. We clarify how and why this is possible 
once the dynamics of the colored Polyakov loop is taken into account. Also we succeed in giving an 
analytic expression for Tc for the transition from two-flavor pairing (2SC) to unpaired quark matter 
in the presence of the Polyakov loop. 
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I. INTRODUCTION 



Quantum chromodynamics (QCD) is expected to exhibit a variety of phases depending on the temperature and 
[ on the baryon density At zero density and finite temperature the two main features of QCD arc the confinc- 
I mcnt/dcconfinement phase transition and chiral symmetry restoration. They should be realized when the hadronic 
Qv^ [ system is heated, for example in ultrarelativistic heavy ion scattering processes such as at the Relativistic Heavy Ion 
. Collision (RHIC) experiment 0] or in the future ALICE experiment at the Large Hadron CoUider (LHC) at CERN. 
Cn| ' Moreover this behavior is clearly seen by lattice QCD simulations [1]. In these conditions quarks and gluons should 
, be released as active degrees of freedom at some critical temperature. Moreover in the same range of temperatures 
' one expects the restoration of chiral symmetry, whose spontaneous breakdown is known to play a key role in the mass 
00 ', spectroscopy of zero density QCD 0, Q. At finite baryon densities similar transitions are also expected although 
' the comparison with lattice simulation data is still lacking due to the so called fermion sign problem. However, 
^ ' at extremely high density, where perturbative techniques are allowed, it is now theoretically well established that 
quarks are deconfined forming diquark condensates so that the system is in a color superconducting ground state with 
asymptotic color-flavor locking (CFL) @. While difficult to achieve in the laboratory color superconductivity might 
be relevant to the inner structure of compact stellar objects Q. 

Exploring phase structure at intermediate density, where neither lattice simulations nor perturbative calculations 
can be trusted, is the object of various model studies. There are different effective models which provide simple 
descriptions of chiral symmetry restoration at finite temperature and density; the Nambu-Jona Lasinio (NJL) model 
is one of them [1, 0, H, Q. The NJL model realizes the spontaneous chiral symmetry breaking of QCD at small 
temperature and density. Despite its simple structure, it can also realize a CFL phase at the largest density. Moreover 
it can even reproduce the correct ratio of the gap and critical temperature for the transition from the CFL to the 
unpaired phase. 

The main defect of the NJL model is the absence of the confincment/deconfinement transition. A theoretical 
attempt to understand the nature of the deconfinement transition goes back to the work [l0| in which dcconfincmcnt 
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was shown to be associated with the spontaneous breaking of the global Z(iVc)-symnietry of a finite temperature 
SU{Nc) pure gauge theory. The order parameter is the traced Polyakov loop, whose condensation and correlation are 
related to the free energy of static quark and the string tension between two static quarks in a thermal medium. 



The inclusion of the Polyakov loop dynamics into the NJL model was first done by Fukushima 11| in order to study 
the relation between chiral restoration and dcconfincment. It is now called "Polyakov loop extended Nambu-Jona 
Lasinio" (PNJL) model. In this model, the chiral condensate qq serves as an order parameter for the chiral transition, 
while the traced Polyakov loop $ performs this job for the dcconfincment transition. Even though the former and the 
latter have their definite meanings, as order parameters, only within different limits, (niq and oo), they 

are still useful as indicators of both crossovers and/or transitions. In addition, the model enables to interpret nicely 
some bulk properties of matter observed on the lattice on the field theoretical ground [l^l . 

The purpose of this work is to investigate the color superconducting phase structure in (T, mg//i)-plane within the 
PNJL model, and to study the effect of imposing neutralities, in both the paired and unpaired phases, in presence of 
the Polyakov loop. The neutrality constraints are known to be important for the candidates of color superconducting 
phases at a realistic density; they open a window in a phase diagram to intriguing gapless phases [13, 13| • Although 
a few works have already explored the pairing phases of PNJL models [l^, [ly, [l^ , none of them takes into account 
either the possibility of complicated gap structures or the neutrality effects. Thus our study is a natural extension of 
them. One surprising result is that once PNJL model is applied to finite density, it inevitably lacks color neutrality 
even when the system is unpaired. This is a sort of sign problem at finite density. In this paper we still proceed on 
the heuristic assumption that PNJL is applicable to finite densities. Also it will turn out that the inclusion of the 
Polyakov loop greatly affects the phase diagram by stabilizing the two-flavor pairing (2SC) phase, and it also brings 
about a color-flavor unlocking transition |18| at finite temperature in a new mechanism. 

The paper is organized as follows. In Sec. |TT1 we introduce our model and approximations. In the first part of 
Sec. mil we demonstrate the lack of color neutrality in the conventional PNJL model at finite density. The rest of the 
section is devoted to discussion of the numerical results. The neutrality effects, the effect of dynamics of Polyakov 
loop, and their interplay will be particularly covered. We summarize the main contents of our paper with some 
concluding remarks in Sec. IIVI 

II. FORMALISM 

In order to accommodate for pairing in the — 0+ channel at finite density, we add the 4-point vertex to the free 
part of the Polyakov-quark model, which hereafter we shall refer to as the Polyakov NJL (PNJL) model. 



Cash 9; ^4] = + 7o(/i + ^Moff))9 + jqP^fq'^Pr,q - U{T, $[^4], $[^4]*). (1) 



Here q stands for the quark field, and a summation over color and flavor degrees of freedom is understood. Prj = 
Cj5erjijerjab = To^iJto) IS the matrix, antisymmetric in color, flavor and spin, specifying the pairing channel. 
The constant G parameterizes the strength of the coupling leading to diquark condensation. We work within the 
chiral SU{2) limit, setting m„ = rud = 0, and take into account the strange quark mass within the high de nsity 
approximation. This means that we include the effect of its finite value in the chemical potential difference Sn^ff UM- 
As a result, 
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S^J.cS = -^J■eQ + Ms^s + MsTs - ^diag.(0, 0, 1)/ x 1^, (2) 

where Q = diag.(2/3, — 1/3, — 1/3)/ x 1^, T3 = 1/ x ^Aa, and Tg = 1/ x ^^Ag, with {Aq} being the standard 
Gell-Mann matrices. We find it more transparent switching to a new spinor basis for the quark field defined as 

QA = {qur,qdg,qsb, QugiQclr, qsr,qub, Qdb, Qsg) by mCauS of: 

9 

qia = '^{FA)iaqA (3) 

A=l 

with {F^jia unitary matrices in color and fiavor space defined in [Toj . In this new basis ^ takes the form of a 
diagonal matrix: J/x^Jab with A = l(wr), 2{dg), • • • , 9{sg). 

We treat the Polyakov loop by the static, homogeneous, and classical background gauge field A4 = igAQ^ where 
the temporal gauge field A4 is introduced by parameterizing the Wilson-line, as L = e'^-^/^. In the PNJL model 
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one assumes that this background gauge field couples to quarks with covariant derivative P„, = — S^qA^. In the 
convenient gauge called Polyakov gauge, the Wilson line L is in the diagonal representation i.e., 

Moreover we restrict ourselves to the case (f>s = such that the traced Polyakov loop $ = ticL/Nc becomes real 
[TH whereas at finite density there is no strict reason why $ should be real [lO, [HI • Thus in this representation, 
(f> — ^'=°'^('^3/T)+i ^ ^j^g effect of the background field A4 is just to shift the color chemical potential to the imaginary 
direction fJ-3 ^ fJ-3 — 2i(f>3 = jl^. For the Polyakov loop effective potential hi we use the following form, inspired by the 
strong coupling analysis of the pure gauge sector [2^, 



rp4 

with 



$*$ + b{T) log (1 - 6$*$ + 4($*-^ + - ($*<I>)^) , (5) 



62(r) = ao + ai(^^)+a2(^^y, 6(r) = 63(f)'. (6) 

The numerical values for coefficients are determined by fitting several quantities to the lattice results of pure gauge 
theory [i2|: 

ao = 3.51, ai = -2.47, 02 = 15.2, ,63 = -1.75. (7) 

In the absence of dynamical quarks, Tq is set to the value of the transition temperature for deconfinement, i.e., 
To = 270 MeV. In our model, we use the value To = 208 MeV which is the theoretically suggested value for To in the 
presence of two light flavors, Nf = 2 [l^, [23|, although as we treat the strange quark mass as a free parameter, our 
calculation will cover the situations between two flavor and three flavor, i.e., Nf = 2 (+1) — > 3. In the case Nf = 3, 
the slightly lowered value To = 178 MeV is proposed, but we have checked that choosing this value for Tq does not 
change our results in any signiflcant way. 

By introducing a charge conjugated fleld qc = —Cq* as an independent field, and after introducing the Hubbard- 
Stratonovich field A^{t, x) = ^Q^^vl ^"^^ ^r]{T, x) = ^-qP^q^ , we integrate out the fermion field Q = {q, qc). Within 
the mean field approximation for A and A, the effective potential becomes 

n(A,„^4,/^e,/i3,M8;/i,T) = Z^(T,<i>,$*)-^-^-Z£!ll 

where £^^ [^4] = I?^[^4]|ao^ia^ with r denoting the imaginary time, and the transpose operation X*" only acts on the 
color and flavor structure of X . 

It is useful to write down the thermodynamic potential in the A,, = and /ie.3,8 = case. Leaving the trace over 
color, it takes a form 

n = U{T, $*) - 2iV/T / -^tr.ln \{l + L'fe~(''-^)/^)(l + Le-(f+^)/^)l . (9) 

J {2ttY L J 

with Nf = 3. Within the Polyakov gauge and imposing the (ps = prescription, the Wilson line takes the following 
form 

I* I , (10) 



^ One can always find the gauge rotation U such that ULU~^ becomes diagonal. This is simply a gauge fixing, and any physical quantities 
will not depend on the gauge freedom U so we can safely reduce the eight dynamical variables to parametrize L up to two independent 
parameters {(/)3,<7ig}. We note, however, when the diquark condensation is taken into account, this is no longer justified unless more 
general ansatz for the diquark condensation, A^^/e^^je^/^j, is adopted. Simultaneous color-flavor rotation can make the condensate 
matrix diagonal such that the usual assumption, A^^/ oc S^^i ,is recovered, but this is nothing but the gauge fixing. Thus in principle, 
if we adopt the diagonal ansatz for the diquark condensation, we can no longer make L gauge rotated to diagonal, and on the other 
hand, if we adopt the diagonal form of L, we should work in more general assumption for A^^/ . Nevertheless, we work in the simplified 
assumption that both L and A,|^/ are of diagonal as in [Tsll leaving a further detailed analysis in the future. 
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with I = e'-^^/^. 

Let us now consider the case with 7^ and fie.a.s 7^ 0- We try to simphfy the expression for the thermodynamic 
potential. Within the current approximation (treating 7^ as a shift of the chemical potential), the action does 
not mix the left-handed quark and the right handed quarks. Thus, we can rewrite the functional determinant as 

^K-^Jy,JnomA,]+^^ + 6f.^^ -A,e,,,e , \ _ ^ a a,). (11) 

Now the Dirac gamma matrices can be regarded as two dimensional matrices for Weyl spinors, say, 7^ = (1,<t). 
Finally, putting p = fiv + 1 with velocity \v\ = 1, and discarding the antiquark contributions, we get the high density 
effective theory (HDET) approximation for the effective potential [l^ . 



where lJc is a momentum cutoff, and the positive energy left handed projected propagator is defined as 



(12) 



^ ' \ -A^e^abe^y iuj„ + v-l-Sfilff + iAl J ^ ' 

This is now a 18 x 18 matrix defined in the color-flavor space. It has a form (iWnlig — Ti). In order to evaluate the 
Matsubara summation, we have to evaluate all the eigenvalues of the hamiltonian density Ti.. Since we have doubled 
the degrees of freedom by introducing the Nambu-Gorkov notation, the eigenvalues of the hamiltonian will appear 
as {Ea{1\\),—Ea{1\\)} with A = 1, 2, • • • , 9. In contrast to the standard NJL models without Polyakov loop, Ti is 
no longer Hermitian due to the imaginary chemical potential jl^, and accordingly each quasiparticle energy Ea{1\\) 
can take in general complex values. Consequently, f2 is no longer restricted to be real. We avoid this sign problem 
by taking the real part of fl as in [l2|. Once these quasiparticle energies are evaluated, and choosing the basis as 
{Ea,—Ea} such that both the conditions, ^Ea > and Ea — i^M^Jsl when A^,(/)3 — ^ 0, arc satisfied, we can 

perform the Matsubara summation as 



miA„ ^4, Me, m, Ms; M, T) = UiT, cO, $*) _ ^ _ _ Z^L _ J^^^^ 



fi" dl II 



+ ELi ^ ["fiEAih) - \h - ^Me1f I + 2Tln(||l + e-^-('ii)/^||) 

(14) 

We took the energy density in the vacuum without any condensation as the reference energy density. The fifth term 
is nothing but the zero temperature part of free quark contribution to the effective potential. We need the ultra violet 
cutoff uJc only in the integral representing the condensation energy, i.e., the last term.^ Instead of G, we use Aq, the 
zero temperature CFL gap in the chiral SU (3) limit without Polyakov loop, as the indicator of diquark attraction [2^ 



G 7r2 \2^/3Ao 

With the use of this cutoff dependent couphng constant, the derivatives of the effective potentials, g^^^" 3) , now have 
well-defined limits as ujc 00. In this way, we can remove the cut-off dependence from the gap equations, while it 
remains in the effective potential itself. 

The evaluation of the effective potential is carried out by finding the eigenvalues of Ti. for given momentum and 
integrating them over the momentum. Then the mean field solution for the ground state is obtained by minimizing 
the effective potential with respect to (A,,, ^3) imposing the proper constraints of charge neutrality, g(^^^^^g) = 0. 



III. RESULTS 

In this section, for the numerical computations we fix = 500 MeV, and concentrate on the case with Aq = 60 MeV 

2 

with LOc ~ 300 MeV. We would study the phase diagram, and how the physical quantities behave as functions of -2^, 



^ The integral of the thermal part {Tln(- ■ • )) can be evaluated without cutoff. 
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and temperature T, treated as free parameters. In the numerical calculations, in order to take the minimum strong 

coupling effect into account, we bring back the momentum dependence of the density of state -^-^^^r^- By 

doing this, the particle-hole asymmetry which is known as the first correction to the weak coupling approximation 
will be properly taken into account. 



A. Color neutrality 

Before discussing in detail the calculation let us briefly comment on a strange feature that shows up, even in 
the unpaired phase (A,; = 0), if one assumes vanishing charge chemical potentials. The interesting fact is that the 
conventional PNJL model calculation at finite /x lacks color neutrality. It is worth stressing here that this does not 
follow from HDET approximation (see ([9]) obtained before introducing HDET formalism). It is clear from (jlOp that 
the Wilson line couples to each color of quarks, (r, g, b), with different weight. Looking at the real part only, the weight 
for (r, quarks differs from that of b quarks. Thus, the energy required to populate {r,g) quarks in the background 
gauge field A4 is different from that for b quarks. Furthermore, cos(03/r) can take a negative value close to —1/2 
when the system is nearly confined. This means that thermal excitations of on-shell (r, g) quarks reduce the pressure 
of the system making it difficult to create on-shell (r, g) quarks in the system. This can be regarded as the effect of 
confinement in this PNJL model at finite density. The problem is that, even in this situation, b quarks can be excited 
at finite temperature because L33 is unity and thus the thermal weight for b quarks does not differ from that in the 
deconfined phase. This unphysical feature might be viewed as an artifact of this PNJL model which originates in the 
assumption that the dynamics of traced Polyakov loop $ can be equivalently described by the constant background 
gauge field A4. As a consequence, the development of the finite value of A4 breaks not only the Z(3) center symmetry 
but also the color SU{3) symmetry in a spontaneous way; this fact does not contradict the Elitzur's theorem [23| 
since we still expect physical quantities such as quasi-particle dispersions should not depend on the gauge and thus 
have their definite meanings even after gauge unfixing which would make A4 itself vanish. 

One may think that this undesired feature is just due to the wrong ansatz of the PNJL model itself. Apart from 
such a possibility, in order to avoid an unphysical appearance of color density within the model, we should inevitably 
introduce an appropriate color chemical potential. According to the gauge of the Polyakov loop, in our case, /is is 
required to maintain the unpaired phase color neutral. To be more explicit, ignoring the antiquarks we can write the 
color density as 



(16) 



nr,g = 2Nf J j^fpip- - {I ^ -^i) 

I OAT r dp 3($-l)tanh(£^) 

where ni, = 2Nf J ^^^y^ ^(p-p^)/t_^i is just the density of a free Fermi gas and 

MP-f^-J)^n{ ^_,^,^_ly^^^ ) (17) 

is the modified Fermi distribution in the presence of A4 describing (r, g) quarks which is plotted in Fig.[IJb) compared 
to the standard Fermi distribution frip— jj) for blue quarks. From it follows that the density of r and g quarks 
differs from that of simple fermi gas, and the difference never disappears unless either T— >Oor<I>^l(r— > 00) is 
approached. The difference also cancels at /x = thanks to the equal and opposite contribution from antiquarks. 

In Fig. [ija), we plot the ratios of each color density, n^, rig, and rif,, to the total quark density as a function of T in 
the unpaired phase with /ie,3.8 — 0- In the numerical evaluation we set /i = 500 MeV, and ignore the tiny antiquark 
contributions. As discussed above, we can see the finite difference between = Ug and njj in the intermediate range 
of temperature."^ The deviation becomes maximum near the steepest point of the Polyakov loop, which is usually 
identified as the deconfinemcnt transition Il2|. 



^ It should be noted that the color density itself is a gauge dependent quantity and thus should depend on the choice of the gauge 
parametrizing the Wilson line. With our diagonal representation of Eq. ^ with ijig = 0, the Tg color density becomes finite as we 
observed above. If we selected the different gauge, the other entry of octet color densities {{q^Tc,q)} should have appeared. The 
important thing is, however, whichever gauge we choose, some color density should become finite; in fact the squared sum of the octet 
color densities is shown to be the gauge independent quantity [2^. 
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FIG. 1: (a): Each color density of a finite density PNJL model with rUs = — /ie,3,8 = as a function of T. Tg-color density 
is induced by the antitriplet color charge I* of Wilson line coupled to quarks, (b): The occupation number profile of quark with 
red, green, and blue, as a function of momentum at T = 0.4ro which is indicated by the vertical dotted line in (a). 



B. Phase diagrams 

In Fig. [2Ia), the phase diagram in the (t, -plane of the charge neutral PNJL model is displayed. The bold 

line represents the first order phase transition, while the thin line represents the second order phase transition. The 
dashed line corresponds to the crossover from the 2SC phase to the gapless 2SC (g2SC) phase. For comparison, we 
have shown in figure[2{b) the phase diagram for the charge neutral NJL model which is the same reported in [2^. We 
are specifying each phase as summarized in TABLE U] From these graphs, the impact of the Polyakov loop dynamics 
on the pairing is quite obvious. The inclusion of the temporal giuon field significantly broadens the region for the 
superconducting phase, in particular for the 2SC phase. In fact, the critical temperature is almost twice as large as 
that in the NJL model without the Polyakov loop which is already known in [l^ . (Note that the scale for the T-axis 
of figure (a) is twice as that of figure (b).) Apart from this significant quantitative change, the qualitative behavior 
of the phase diagram is not so much affected. In both cases, the d-quark superconducting phase (A2 = 0; dSC) exists 
in a small region at finite temperature [1^, HO], and there is the doubly critical point indicated by the upper triangle, 
a point where the line for the vanishing of Ai intersects that for A2 [2^. Also the existence and the location of the 
critical point where the fully gapped CFL phase turns into the gapless CFL (gCFL) phase phase is not affected. In 

both figures the point is indicated by a lower triangle on the ^ axis. This fact means the the effect of the Polyakov 
loop on the pairing is absent at T = 0. The reason is that at T = the Polyakov dynamics decouples from the pairing 
(NJL) dynamics so that the effect is absent because the temporal gauge field (p^ is proportional to T itself. 

C. Impact of the Polyakov loop at finite temperature; color-flavor unlocking and stiffs 2SC phase 

Next we focus in detail on the phase transitions at TOs = in order to study the impact of the Polyakov loop 
dynamics and charge neutrality. To this end, we examine each effect step by step. In Fig. [3l^a), we show the gaps 
A^(r) (solid lines) and the Polyakov loop <i>(T) (long-dashed one) calculated without the neutrality. For comparison, 
we also show by the dashed line the A,,(T) calculated with the NJL model without the Polyakov loop. In this case, the 
three gaps have the same behavior as functions of T and they drop to zero simultaneously when ^ 0.714Ao (shown 



(g)CFL 


A, 7^ for = 1, 2, 3 (with A2<\^ + ^ + ^-^ 


orAi< -ia + i^_^ satisfied) 


dSC 


Ai,3 / 0, A2 = 


uSC 


A2,3 /O, Ai = 


(g)2SC 


A3 / 0, Ai = A2 = (with A3 < 


satisfied) 



TABLE I: The definition of pairing phases of current interest. 
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FIG. 2: (a): Phase diagram in (^,T)-plane at Ao = 60MeV, /i — 500 MeV with the Polyakov loop under charge neutrahty 
constraints. The energy scales are normalized by Aq. The bold line corresponds to 1st order phase transition while the thin 
line indicates the 2nd order phase transition. For other instructions, see text, (b): The same as (a) but without the Polyakov 
loop. In both figures, the bold square put on the T axis indicates the weak coupling approximation of the critical temperature 
at chiral limit without $, i.e., T"/Ao = 0.714. 




FIG. 3; (a): Gaps (solid lines) and the Polyakov loop (long dashed line) as functions of T without the neutrality constraints, 
i.e., Me, 3, 8 = 0. For the dashed line and the bold square, see the text, (b): The same as (a) but with the neutrality constraints 
being respected. Three charge chemical potentials are also depicted. A3 (ltd) in (a) is also shown by the dashed line, just for 
comparison. 



by the bold square on the ^ axis) is approached. Once the Polyakov loop is taken into account this is no longer 
true as one can sec from the figure. Two gaps have the same magnitude, Ai{ds) — A2(.su), while A3{ud) is larger; 
moreover the two gaps Ai ~ A2 drop to zero simultaneously near a point lower than the bold square (T^ ~ 0.714Ao). 
This can be described as a second order color-flavor unlocking transition induced by the Polyakov loop dynamics. 
This behavior is not strange because even though the Polyakov loop is blind to the flavor degrees of freedom, in 
the color-flavor locked phase, however, color is locked to flavor which is the way the Polyakov loop affects the gap 
structure. In fact, as we have already discussed in Sec. IIII A[ the presence of the Polyakov loop induces a finite color 
Tg density as = 7^ rib in the unpaired phase. This means that the existence of the Polyakov loop adds to the 
real part of the effective potential, 3?J7, the finite external field with Tg charge, which explicitly breaks global color 
SU{3)c down to SU{2)c- Since color and flavor are locked in the CFL phase, the external field induced by $ tends to 
break the SU(3)c+v symmetry down to SU{2)c+v- This is the very reason for the splitting, Ai = A2 7^ A3 and also 
for the color-flavor unlocking to the 2SC phase. The emergence of color-flavor unlocking is one of the most interesting 
features of the inclusion of the Polyakov loop dynamics in NJL model. 

A further aspect deserves to be stressed. The 2SC phase persists up to I.SAq which is almost twice as large as 
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the weak coupling formula for the critical temperature, T° = 0.714Ao (indicated by the bold square). This striking 
feature is already noticed in the limiting case with the pure CFL ansatz using the Ginzburg-Landau approach p^ . 
In what follows, we try to explain this fact. In the case of the 2SC case, we can explicitly determine the quasiparticle 
energy dispersion 



Four out of nine quasiquarks have nontrivial dispersion laws: -©1,2 ~ Ti</ 



A3 and 



i?4.5 = ±iCf 



A|. In this case, there is no shift of the averaged chemical potential to the imaginary direction. 



The gap equation in the absence of charge chemical potentials then becomes 



din 



1 - 2ff 



1 - 



|(l-<E.)(l-2/f)(l-/f.) 



1-3(1-<E.)/j=.(1-/f) 



(18) 



is the Fermi distribution, and its energy argument is now supposed to be ^/l 



A|. The second 



where frix) = 

term proportional to fp is referred to as a blocking integral due to thermally excited quarks. We need no cutoff ojc 
for this part to be evaluated. We notice that the effect of the temporal gauge field is just to suppress the blocking 
integral. This explains the robustness of 2SC in the presence of (^3, and the increase of the critical temperature as 
follows. We put A3 = in (fTSl) and try to solve it in T to derive the critical temperature. (US]) with A3 = is nothing 
but the condition for criticality, the Thouless criterion, which guarantees the divergent susceptibility. This condition 
together with the definition of the effective coupling coupling constant ([TS]) leads to 



In 



T 



eiB 21/3 Ar 



din — tanh 



k_\ 3(1 - $)/f (1 - fp) 
2t) l-i{l-^)fF{l- fp) 



^(1-$). 



(19) 



Note that the quantity J- is dimensionless and does not depend on T. In the case of the deconfinement phase with 
$ = 1, JF vanishes so that it simply reproduces the standard expression for the critical temperature T° = 2^2^/"^ Aq = 
0.714Ao. It is now easy to imagine that the deviation of <& from unity leads to the positive and thus increases the 
critical temperature = 0.714Aoe-^(^^*^. To the first order in (1 — $), can be calculated as 



T{1 - $) 



21C(3) 

47r2 



1 - $) 



0.64(1 - $) 



Thus to this order, the critical temperature is approximated by 

= 0.714AoX e"-64(i-*), 



(20) 



(21) 



near <& 1. When $ = 0.4 is substituted into the above formula, gets the factor of enhancement e^-^^i^-"^) ^ 1.5. 
Although it is within the linear level, this value fairly agrees with the numerically obtained factor, 1.8. If we use the 
numerical value of T{\ — 0.4) = 0.584, the factor of enhancement becomes 1.79; the agreement is perfect. Although 
unrealistic, at $ = (the confinement), the analytical evaluation is also possible. In this case we have = InS'^/^ 
so that we have the factor e-^ = 3-\/3 = 5.2. This is the theoretical maximum of the critical temperature in the PNJL 
model at weak coupling. 



D. Effect of charge neutrality; the two-step hierarchical unlocking transition 

Let us now discuss charge neutrality at TOs = 0. In Fig. [3Kb), we show the gaps and chemical potentials as a function 
of T calculated respecting the charge neutrality constraints. For comparison, we have shown A3 without neutrality 
(in Fig. [3lja)) by a dashed line. At a first glance, we notice that, even quantitatively, the charge neutrality plays only 
a minor role at mg = 0. 

However, several interesting remarks deserve a discussion here, (i) First, the charge neutrality conditions lift the 
degeneracy Ai = A2 away and open a small window for the dSC (A2 — 0) phase between the CFL and 2SC phases, 
(ii) Second, ^3 vanishes when the 2SC phase sets in. (iii) Lastly, //g does not vanish even when all the pairing melt 
and the system goes into the unpaired quark matter as already discussed in Sec. IIII AI in the case of no pairing at 
all. Without finite /xg, the unpaired system inevitably has a finite Tg-color charge. The reason for (ii) is simple. The 
2SC pairing preserves the SU{2)c symmetry intact or, in other words, the 2SC gap has SU{2)c singlet structure and 
is transparent to the SU{2)c charge. Therefore the system with no Ts-color charge should have ^3 = 0. One may 
think that a finite value of 03 induces a finite T3 charge in the system with /Z3 = 0, but this is not correct. In fact, 
as we saw in the discussion in Sec. IIII A( restricting ourselves to the real part of fl, Ur = rig ^ rib is realized in the 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

(m//2ii)/A„ (m//2ii)/A„ 



FIG. 4: (a): The gaps A,, as a function of at T = 0. The dashed hne is A3 for the 2SC solution. The vertical dotted 

line corresponds to the 1st order gCFL-to-g2SC transition; on the left of it the CFL is realized, while on the right of it the 2SC 
is realized as the ground state, (b): The same as (a) at finite temperature, T = 0.3Ao. All the transitions are of second order. 

unpaired system inducing a Tg charge but no T3 charge. The point (i) can be understood as follows. In the presence 
of 03, the gaps split as Ai = A2 < A3 as we saw above. This induces an imbalance in the thermal population of 
nine quasiquarks in the CFL phase. Accordingly, the charge neutrality is lost unless /is, /ig and are tuned to their 
appropriate values. But of course the finite values of and /is explicitly break the remaining SU (2)c+y symmetry. As 
a result of this secondary effect, the SU{2)c+v degeneracy should be lifted away, as Ai 7^ A2. The appearance of the 
dSC phase at rus = is in contrast either to the Ginzburg-Landau approach [30| or to the pure N JL calculations [2^ ; 
this is definitely due to the nontrivial interplay between the neutrality constraints and the Polyakov loop dynamics. 

E. Effect of the stress due to nonzero strange quark mass 

Let us now discuss the effect of a nonzero strange quark mass on the structure of gaps. In Fig. HJa) we show the 

2 

zero temperature gaps A,, as functions of both for the CFL (solid lines) and 2SC (dashed line) solutions. The free 

energy comparison shows that there is a first order phase transition from the gCFL phase to the g2SC phase at the 

2 

point ^ = 2.4Ao indicated by the vertical dash-dotted line in the figure. Since at T = the Polyakov loop dynamics 
decouples from the pairing (NJL) sector, the Polyakov loop plays no role in the gap structure. Consequently, the 
phase structure and the behavior of the gaps is similar to the result of [2^ although the HDET approximation was 
not adopted there. 

At low rris, the CFL phase is realized, and it continuously goes into the gCFL phase [l^ at a point slightly lower 
2 2 
than ^ ~ Aq. Then eventually the gCFL phase is taken over by the g2SC phase at ^ = 2.4Ao by a first order 

transition [13, H^. In Fig. [4Ub), we have shown the gaps A,,(^) at a finite temperature, T = O.SAq. The first order 
transition is completely washed away, and there are two successive second order unlocking transition until the system 
gets unpaired, first from the CFL phase to the uSC phase, and subsequently from the uSC phase to the 2SC phase. 
This feature is also qualitatively the same as in the NJL calculations [29j . 

F. Interplay of the Polyakov loop dynamics and enforced neutrality at finite strange quark mass 

Let us finally examine the impact of enforcing charge neutralities and including the Polyakov loop dynamics into 

2 

the pairing phases at an intermediate density represented by the finite value of In Fig. [^a) the gap A3 and 

2 

the charge chemical potentials {/.te,/i8} are depicted as functions of T at ^ = 3.25Ao. At this value of the stress, 
the CFL pairing is not possible, so only the (g)2SC phase can show up as a pairing pattern. What is surprising and 
also intriguing is that A3 once melts at T = O.IAq but appears again at higher temperature about 0.24Ao, and then 
finally vanishes completely when T exceeds 1.23Ao. The 2SC phase exists in two different region in temperature. 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 



T/A„ T/A„ 

FIG. 5: (a): The gaps A,,, the chemical potentials fj.e,s, and the traced Polyakov loop $ as a function of T at — 3.25Ao. 

(b): Comparison of (a) with simplified versions of the model. The bold full line is the same as A3(r) in (a). The long-dashed 
line is the 2SC gap A3(T) calculated in the pure NJL model without the Polyakov loop dynamics. The dashed line is A3(T) 
calculated with the Polyakov loop but without respecting charge neutralities, i.e., ^e,3,8 = 0. 



This feature is definitely due to the inclusion of both the neutrality and the Polyakov loop into the problem. To see 
this, we show in Fig. [SJb), the gaps calculated with simplified versions of our PNJL model. The solid line indicated 
by "Full" is the same as A3 in Fig. [5l[a). The long-dashed line represents the result calculated using the pure NJL 
model without the Polyakov loop dynamics, while the dashed line is that calculated with PNJL model but without 
imposing the charge neutrality constraints, i.e., putting /ie.3,8 = from the very beginning. From these comparisons, 
it is clear that the appearance of the intriguing possibility of the existence of two islands of 2SC in temperature is 
due to the combinatory, cooperative effect between the Polyakov loop dynamics and the neutrality constraints. In 
contrast to the case with mg = 0, imposing neutrality has a sizable effect on the gap. It significantly reduces the 
magnitude of the gap. It is so because in the case of a finite value of m^, not only /xg but also fie should be finite even 
in the unpaired phase in order to guarantee electrical neutrality. Moreover the effect of the Polyakov loop is not only 
to stabilize the 2SC phase against the increase of temperature as in = case, but also to suppress the pairing at 
low temperature making two separated islands of 2SC in temperature. 



IV. SUMMARY 

In this paper, we have studied the quark matter phase structure in [t, plane starting from the PNJL model 

in which a temporal static gluon field couples with quarks. This work is a natural extension of previous studies 
[isl . ITgI . flTj . The particular focus was put on the effects produced by the inclusion of the Polyakov loop dynamics on 
the pairing phases and by the enforcement of color and electrical neutrality. 

In the conventional PNJL model, there is a mismatch in each color density so that the model lacks the color 
neutrality even in the unpaired phase. This unphysical feature is significant in the proximity of the deconfinement 
transition. We have pointed out that this behavior may be due to the original assumption hidden in the PNJL 
model i.e. that the traced Polyakov loop dynamics can be represented by the inclusion of the static temporal "colored" 
gauge field which couples to the fundamental color charge of dynamical quarks. By this assumption one misses gauge 
invariance. Once this fact is admitted, in order to avoid the unphysical appearance of color densities within this 
model, one has to include the charge chemical potentials into the problem from the beginning. In fact, we have shown 
that /ig should be finite to maintain color neutrality in the unpaired phase. 

In the detailed numerical analysis, we have depicted the phase diagram in (T, m^/2/x)-plane, and clarified how the 
phase diagram is affected by the inclusion of the Polyakov loop and the enforcement of charge neutrality. Even at 
iJis = 0, the effect of the Polyakov loop is remarkable; it breaks the SU{3)c+v down to SU{2)c+v and causes a 
continuous color-fiavor unlocking at finite temperature in a novel mechanism. In addition, it makes the 2SC phase 
much more robust against the increase of temperature. The critical temperature is about twice as large as the weak 
coupling prediction, which is consistent with previous calculations [l^, [H, [l3|- We have also examined a formal 
explanation about these facts and derived an analytical expression of this enhancement factor; it turned out that the 
temporal gauge field reduces a blocking integral. 
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The effect of imposing neutralities gives only a tiny effect at rus = although it opens a small window for the 
dSC realization between the 2SC and CFL phases by lifting away the SU {2)c+v degeneracy. In this case, we have a 
hierarchical unlocking, CFL dSC~> 2SC, until it eventually goes into the unpaired phase. This is contrast to the 
pure NJL calculation without the Polyakov loop where the dSC never shows up at = [1^. 

The sizable effect of imposing charge neutrality on the pairing phases manifests itself at finite . We have shown 
that the nontrivial, complicated interplay between the charge neutrality constraints and the Polyakov loop dynamics 
at rris ^ produces a thermal rcentrancc phenomenon, as two isolated windows for the 2SC pairing can show up on 
the temperature axis. 

There are several ways to extend our current study. One is to take the chiral condensation into account by including 
the chiral condensate and removing the high density approximation (3l] |. By this improvement, one can study the 
interplay between the chiral condensate, the Polyakov loop, and color superconductivity at the same time. The 
other possibility is to study mesonic modes [s^ . [ssj as well as the Meissner masses in the superconducting phases. 
This might have an impact either on a possible meson condensation in superconducting phases or on the so called 
chromomagnetic instability problem in gapless phases [H, Is^l • These studies may be presented elsewhere in future. 

To conclude let us stress that two alternatives have presented to us: 

a) application of the PNJL model to finite density is pathological and should be avoided, in relation to the fact 
that color neutrality is not visibly satisfied; 

b) the model can be used also at finite density provided neutrality is enforced: for such a case we have derived the 
detailed consequences obtaining surprising results but without apparent physical inconsistencies. 

We hope that further work will illuminate on the choice between a) and b). 

We thank Kcnji Fukushima, Mci Huang, Andreas Schmitt and in particular, Igor Shovkovy for enlightening discussion. 
One of us (H. A.) thanks I. N.F.N, for financial support. Numerical calculations were carried out on Altix3700 at YITP 
in Kyoto University, and on the workstation NETCLUS at University of Bari (I.N.F.N., Sczione di Bari). 
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